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TITLE OF THE INVENTION 
METHOD AND DEVICE FOR BACKSCATTER SPECTROSCOPY 

CROSS-REFERENCE TO RELATED APPLICATIONS 
5 [0001] This application is a 37 1 of international application number 

■ 

PCT/NL2004/000657, filed on September 22, 2004. 

BACKGROUND OF THE INVENTION 
1. FIELD OF THE INVENTION 

1 0 [0002] The present invention relates to a method of determining of a physical feature of a 

medium, comprising: producing radiation with a light source; placing a probe on a sample of the 
medium, the probe comprising a first optical fiber having a first diameter, and at least a second 
optical fiber having a second diameter; sending light coming from the light source, through the 
first optical fiber; collecting first backscattered radiation through the first optical fiber and 

1 5 second backscattered radiation through the second optical fiber; producing a first signal based on 
the first backscattered radiation, and a second signal based on the second backscattered radiation; 
and determining a measured differential backscatter signal as a function of wavelength using the 
first and second signals. 

20 2. DESCRIPTION OF THE RELATED ART 

[0003] Such a method is known, as described in A. Amelink, M.P.L. Bard, J. A. Burgers, 

and H.J.C.M. Sterenborg, "Single Scattering Spectroscopy for the Endoscopic Analysis of 
Particle Size in Superficial Layers of Turbid Media", Applied Optics 42, pp. 4095-4101 (2003); 
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which describes a special device used to determine particle sizes in superficial layers. The 
device is suitable for measuring particle sizes in, for example, an aqueous suspension with 
polystyrene spheres, but is not fitted to accurately measure particle sizes in living tissue. So, 
determining whether living tissue is normal or precancerous, by way of measuring particle sizes 

5 in living tissue is not very promising. 

[0004] As described in R.M.P. Doornbos, R. Lang, M.C. Aalders, F.W. Cross, and 

H.J.C.M. Sterenborg, "The Determination of In Vivo Human Tissue Optical Properties and 
Absolute Chromophore Concentrations Using Spatially Resolved Steady- State Diffuse 
Reflectance Spectroscopy", Phys. Med. Biol. 44 (1999), pp. 967-981; the optical properties of 

10 human tissue are determined in vivo using a spectroscopic arrangement with ten optical fibers. 
One of the fibers is used to irradiate a sample, and nine other fibers collect the reflected light. 
By using a multitude of fibers to collect the reflected light, it is possible to calculate scattering 
and absorption coefficients of the sample. However, the method is not suitable for locally 
measuring the optical properties of the tissue. In particular, only mean values of the absorption 

1 5 coefficient of a relatively large part of the sample can be determined. 

BRIEF SUMMARY OF THE INVENTION 
[0005] It is an object of the present invention to locally measure a physical feature, such 

as a concentration, of a substance in a medium. 
20 [0006] The object is achieved by a method as described above, characterized by 

calculating the physical feature by curve fitting the measured differential backscatter signal to a 
backscatter function, in which the backscatter function is a function of an average path-length 
traveled by detected scattered photons, wherein the average path-length is independent from an 
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absorption coefficient of the medium, and from a scattering coefficient of the medium. Contrary 
to methods using diffusely scattered photons such as described in Doornbos et al. cited herein, in 
the method according to the invention, the local absorption coefficient of the sample is measured 
in an absolute way, independent of the magnitude of the local scattering and absorption 
5 coefficients. This facilitates the measurement of absolute concentrations of absorbing molecules 
in a sample without requiring prior knowledge of the magnitude of the scattering and absorption 
coefficients of the medium. 

[0007] In an example embodiment, the average path-length is proportional to the first 

fiber diameter. This has as additional advantage that the average path-length and thereby the 

10 average penetration depth into the sample of the photons that contribute to the differential 

backscatter signal can be controlled by choosing the fiber diameter. As a result, the sampling 
volume can be controlled by adjusting the fiber diameter. Hence, the fiber optic probe can be 
engineered to match the relevant dimensions of the medium under investigation. 
[0008] In a particular embodiment, the physical feature is a concentration of at least one 

1 5 substance in the medium. 

[0009] The invention also relates to a device for determining a physical feature of a 

medium, comprising: a light source for producing radiation; a probe with at least a first and a 
second optical fiber, the first optical fiber having a first diameter and being arranged to deliver 
the radiation on a sample of the medium and to collect first backscattered radiation from the 

20 sample, the second optical fiber having a second diameter and being arranged to collect second 
backscattered radiation, wherein the second optical fiber is positioned alongside the first optical 
fiber; a spectrometer for producing a first signal based on the first backscattered radiation, and 
for producing a second signal based on the second backscattered radiation; a processor arranged 
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to determine a measured differential backscatter signal as a function of wavelength using the first 
and second signals, characterized in that the processor is arranged to calculate the physical 
feature by curve fitting the measured differential backscatter signal to a backscatter function, in 
which the backscatter function is a function of an average path-length traveled by detected 
5 scattered photons, the average path-length being independent from an absorption coefficient of 
the medium, and from a scattering coefficient of the medium. 

[0010] Furthermore, the invention relates to a computer program and a data carrier 

provided with a product including the computer program and enabled to calculate a physical 
feature by curve fitting a measured differential backscatter signal to a backscatter function using 

10 an average path-length as described herein. 

[0011] In another aspect of the invention, the invention relates to a method of 

determining a physical feature of a medium, comprising: producing radiation with a light 
source; placing a probe on a sample of the medium, the probe comprising a first optical fiber 
having a first diameter, and at least a second optical fiber having a second diameter; sending 

15 light coming from the light source, through the first optical fiber; collecting first backscattered 
radiation through the first optical fiber and second backscattered radiation through the second 
optical fiber; producing a first signal based on the first backscattered radiation, and a second 
signal based on the second backscattered radiation; and determining a measured differential 
backscatter signal as a function of wavelength using the first and second signals, characterized 

20 by calculating the physical feature by curve fitting the measured differential backscatter signal to 
a backscatter function, in which the backscatter function is a function of a mean free path of 
photons. In this method, it is assumed that only singly scattered photons contribute to the 
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differential backscatter signal and as a result the backscatter function can be easily derived 
analytically. 

[0012] In an example embodiment, the physical feature is a concentration of at least one 

substance in the medium. 
5 [0013] The invention also relates to a device for determining a physical feature of a 

medium with a processor arranged to calculate the physical feature by curve fitting a measured 
differential backscatter signal to a backscatter function using a mean free path as described 
herein. 

[0014] Furthermore, the invention relates to a computer program and a data carrier 

10 provided with a product including the computer program and enabled to calculate a physical 

feature by curve fitting a measured differential backscatter signal to a backscatter function using 
a mean free path as described herein. 

[0015] Finally, the invention relates to a method for simultaneously measuring 

backscatter radiation on different locations of a sample; determining a physical feature for the 
1 5 different locations; and calculating a standard deviation of the physical feature. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 
[0016] The present invention will be described below with reference to exemplary 

embodiments and the accompanying schematic drawings, in which: 
20 [0017] Fig. 1 is a schematic diagram of a measuring device according to a preferred 

embodiment; 

[0018] Figs. 2a and 2b show cross-sections of a sample and two fiber tips in the situation 

wherein the mean free path of the photons is much larger than the diameter of the fibers; 
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[0019] Fig. 3 shows the results of Monte Carlo simulations for a homogeneous medium; 

[0020] Fig. 4 shows a differential backscatter signal normalized at zero absorption for 

several scattering coefficients; 

[0021] Fig. 5 shows a differential backscatter signal of a dilute suspension of 0.2 jam 

5 polystyrene spheres; 

[0022] Fig. 6 shows the total differential backscatter signal as a function of the reflection 

coefficient \i s (\) in the range of 10-100 mm" 1 ; 

[0023] Fig. 7 shows the measured and calculated average path length x as a function of 

the average scattering coefficient; 
10 [0024] Fig. 8 is a graph of measurement results for three different absorption coefficients 

Ha showing the dc-fiber signal I, the c- fiber signal J and the differential backscatter signal R bs as 
a function of wavelength; 

[0025] Fig. 9 shows a typical spectrum of an absorption curve A along with the specific 

absorption coefficient of Evans Blue dye; 
15 [0026] Fig. 10 shows a measured A* as a function of the absorption coefficient jx a at 

X = 600 nm; 

[0027] Fig. 1 1 shows typical spectra measured in a suspension of 1 .0 ^m polystyrene 

spheres with and without Evans Blue dye; 

[0028] Fig. 12 graphically shows a molar extinction coefficient as a function of 

20 wavelength; 

[0029] Fig. 13 shows in vivo measurements and a fit of the differential backscatter signal 

Rb S in a human trachea realized using a fiber diameter of 400 ^im; and 
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[0030] Fig. 14 shows in vivo measurements and a fit of the differential backscatter signal 

Rbs in a human trachea showing very low oxygenation indicative for lung tumor. 

DETAILED DESCRIPTION OF THE INVENTION 
5 [0031] The present invention is related to international application number 

PCT/NL2004/000657, filed on September 22, 2004, which is incorporated herein by reference in 
its entirety. 

[0032] A schematic diagram of a preferred embodiment according to the invention is 

shown in Fig. 1 . The setup consists of a set of optical fibers for the delivery and collection of 

10 light to and from a sample 1 under investigation. Light from a light source 2, for example a 
Tungsten Halogen lamp (Avantes HL-2000-FHSA), is led through a first arm 3 of a bifurcated 
optical fiber. The bifurcated optical fiber is at a distal end 4 coupled to a first distal end of a 
delivery-and-collection fiber 5 (in the following referred to as dc-fiber 5) which is small enough 
to be fit through a working channel of a clinical endoscope, not shown. A second distal end of 

15 the delivery-and-collection fiber 5 contacts the sample 1. Alongside the dc-fiber 5, a collection 
fiber 6 is arranged to collect light reflected by the sample 1 . The collection fiber 6 (referred to as 
c-fiber 6) is connected to a slave channel of a dual-channel spectrometer 7, for example an 
Avantes SD2000. Preferably, the dc-fiber 5 is polished at a small angle to reduce specular 
reflections. 
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[0033] Light reflected back from the sample 1 into the c-fiber 6 is led directly into the 

slave channel of the dual-channel spectrometer 7. A second arm 8 of the bifurcated fiber is 
connected to a master channel of the dual-channel spectrometer 7. Light reflected into the 
dc-fiber 5 is coupled back into the bifurcated fiber, and reaches the dual-channel spectrometer 7 

5 via the second arm 8 of the bifurcated fiber. An output of the spectrometer 7 is connected to an 
input of a processor 9 which is arranged to analyze signals from the spectrometer 7. 
[0034] If only the dc-fiber 5 is used to deliver and collect light to and from the sample 1 , 

a large fraction of collected light is due to single backscattering from small sample depths, as 
described in Amelink et al. cited herein. A single-to-multiple scattering ratio depends on the 

10 scattering coefficient and phase function of the sample 1 and on a diameter of the dc-fiber 5. 
The contribution of multiply scattered light to the signal of the dc-fiber 5 can be approximately 
determined by combining the signal of the dc-fiber 5 with a signal coming from an additional 
fiber, i.e. the c-fiber 6 mentioned above. 

[0035] As described in H.C. van de Hulst, "Light Scattering by Small Particles", Wiley, 

15 New York, 1957; a differential backscatter signal R bs as a function of the wavelength k is 



determined using a formula like 



J'— W) (J^ue M- J black W)y 



0) 



where I(>.) is the signal from the dc-fiber 5 in contact with the sample 1 , l n (K) is the signal from 
the dc-fiber 5 submersed in a fluid with an appropriate refractive index (for tissue: water would 
20 be appropriate), IwhiteOO is the signal from the dc-fiber 5 with the probe-tip at a specific distance 
from a diffuse reflecting reference material with a large, preferably wavelength-independent 
reflectance coefficient (white spectralon) and WiA) is the signal from the dc-fiber 5 with the 
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probe-tip at that same specific distance from a diffuse reflecting reference material with a small, 
preferably wavelength-independent reflectance coefficient (black spectralon). Furthermore, J(X) 
is the signal from the c-fiber 6 in contact with the sample 1 and J w hite(^) and Jbiack(^) are the 
signals from the c-fiber 6 with the probe-tip at the previously mentioned specific distance from 
5 the white or black spectralons, respectively. Finally, c is a calibration constant that depends on 
the distance between the probe-tip and the reference materials. 

[0036] According to the invention, the processor 9 is arranged to calculate the physical 

feature using a predefined mathematical model, the differential backscatter signal (R bs ) and a 
curve fitting mechanism. In an example embodiment, the diameter of the fibers 5, 6 are selected 
1 0 depending on a mean free path (mfp) of photons sent into the sample 1 . It is noted that if the 
mean free path cannot be estimated before selecting a fiber diameter, initially two arbitrary fiber 
diameters may be selected. After curve fitting the measuring results using two different 
mathematical models, it will show which model applies. 

[0037] Figs. 2a and 2b depict fiber tips of the dc-fiber 5 and the c-fiber 6 in the situation 

15 wherein the mean free path (mfp) of photons coming out of the dc-fiber 5 is much larger than a 
diameter fiber of the fibers 5, 6. In an example embodiment, the diameters of both fibers 5, 6 are 
of equal size. However it should be understood that other selections are possible. In Fig. 2a, 
lines 21 and 22 show an example of a path traveled by a detectable singly scattered photon. In 
Fig. 2b, lines 23, 24, 25 and lines 23, 24, 26 show two possible paths of detectable multiply 
20 scattered photons. All multiple scattering events occur at such large distances from the fiber tip 
of the fibers 5, 6 that the probability of detection of multiply scattered photons is roughly equal 
for the dc-fiber 5 and the c-fiber 6. The differential backscatter signal Rb S (X) will now purely be 
determined by singly scattered photons. 



9 



SUBSTITUTE SPECIFICATION 



207,517 



[0038] In an example embodiment, the respective diameters of the fibers 5, 6 are selected 

so that mfp>dfiber. In the predefined mathematical model of this embodiment, the differential 
backscatter signal Rbs(^) is an exponential function of two times the mean free path. Below, an 
explanation for this model is given. 
5 [0039] In the absence of absorbers, the differential backscatter signal Rb S (X) is 

proportional to the local, superficial scattering coefficient M*(^) = Qsca(^)*P*As: 

R b M) = C app ~- ldn p(A,Q) Q sca (A) p As (2) 

"AM 

where C app is an apparatus constant that depends amongst others on the distance between the 
probe tip and the reference materials (black and white spectralon), p(^ft) is a function called the 
10 phase function where Q is the scattering angle, Q sca (X) is the scattering efficiency, p is the 

concentration of substances present in the sample 1 , and As is the area of a scattering particle. 
For example, using a fused silica fiber with numerical aperture NA= 0.22, the differential 
backscatter signal Rb S (X) can be approximated by 

1 In 180 

K W * C app ~-\d<p- \d0 ■ sin(0) • p(A,\ 80) • m s W = C app ■ 80) • ju s (A) (3) 



4^ 0 1 70 



15 where cp is the azimuthal angle and 0 is the polar angle. 

[0040] Fig. 3 shows the differential backscatter signal Rb S (X) of measurements (see dots 

in Fig. 3) of a dilute suspension of 0.2 |am polystyrene spheres along with a calculation (see 
curve 32) according to Eq. (3). In Fig. 3 R bs (^) is shown using arbitrary units (a.u.). Also a 
value of Q ra dar where Q rad ar = 4tc-p(X, 180)- Q sca (>0, is indicated in the figure. Fig. 3 shows 

20 excellent agreement between the measurement (i.e. the dots) and the calculation, which indicates 
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that if mfp>dfiber 5 the single scattering is indeed the dominant contributor to the differential 
backscatter signal Rbs(^) as defined in Eq. (1). 

[0041] A singly scattered photon first travels from the tip of the dc-fiber 5 to a particle, 

and then (the same distance) back from the particle to the tip of the dc-fiber 5 (or tip of c-fiber 6), 
see also Fig. 2a. So an average path length x(k) traveled by the measured single scattered 
photons is equal to two times the mean free path mpf(X), i.e. 

r(A) = 2 • mjp(X) (4) 

[0042] In the presence of n absorbing species with specific absorption coefficients, 

|x a spec,l (A,), the differential backscatter signal becomes 



n 



R bs = C app • pWL 80) • Ms W ■ exp(-r (A) • £ A • nT* W) 
10 w „ (5) 

= C m ■ p(A,l 80) MsW' exp(-2 • mfp{X)) • £ Pi ■ nT' W) 

1=1 

where Capp' is an apparatus constant, p(X,,180) is the phase function, \i s (X) is the scattering 
coefficient of the medium, X is the wavelength of the first and second backscattered radiation, 
mfp(A,) is the mean free path as a function of the wavelength, n is the number of substances in the 
sample 1, pi is the concentration of absorber i present in a detection volume of the sample 1, and 

15 (i a spec,i (^) is the absorption coefficient of absorber i as a function of the wavelength. 
[0043] It is noted that in Eq. (5) the assumption is made that absorbers are 

homogeneously distributed and do not influence each other. The Eq. (5) may be corrected for 
non-linear phenomena such as an inhomogeneous distribution of absorbers, see e.g. R.L.P. van 
Veen, W. Verkruysse, and H.J.C.M. Sterenborg, "Diffuse Reflectance Spectroscopy from 500 to 

20 1060 nm Using Correction for Inhomogeneously Distributed Absorbers", Opt. Lett. 27, pp. 246- 
248 (2002). 
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[0044] According to an example embodiment, the specific absorption coefficients of the 

absorbers, the wavelength dependency of the scattering coefficient |w s and the phase function p, 
together with Eq. (5), are used in order to calculate the concentrations of all the absorbing 
substances present in the detection volume of the sample 1 . Since the detection volume is 
5 typically very small in the present invention, the extracted concentrations are highly spatially 
resolved. This is not possible with the known methods that are based on diffuse reflectance, and 
wherein the obtained concentrations are averages over large sample volumes, see e.g. Doornbos 
et al. cited herein. 

[0045] The apparatus constant C a pp 5 (Eq. 3) can be determined for a specific distance 

10 between the tip of the dc- fiber 5 and the reference materials (black and white spectralon). For a 
suspension of monodisperse polystyrene spheres of known size and concentration, the scattering 
coefficient and the phase functional p(180) can be calculated using Mie theory as described in 
the van de Hulst article cited herein. The apparatus constant C app ' simply follows from Eq. (3). 
In terms of the volume fraction f of the suspension, the radius of the spheres a and the radar 
1 5 efficiency coefficient QradarM = 47t-p(^, 1 80)- Q sca (?0 the apparatus constant is determined by 

K W « C ap ' p • pOU 80) • Ms (A) = C app • 0.05968 ■ L . (A) (6) 

[0046] According to another embodiment, the selected diameter dfiber is chosen so that 

the mean free path is smaller than dfiber- hi this embodiment, the differential backscatter signal 
R bs is a function of the fiber diameter fiber d fib er. This will be discussed in more detail below. 
20 [0047] When the mean free path of the photons is smaller than the selected fiber diameter 

(i.e. mfp(X)< dfiber), the contribution of multiply scattered light to the differential backscatter 
signal Rbs(A,) of the single dc-fiber 5 cannot completely be removed using Eq. (1). In this case, it 
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appears that the average path length of the photons contributing to the signal Rbs(^) becomes 
nearly independent of the optical properties of the sample 1 . In this situation, multiple scattering 
events already occur at small distances from the tip of the dc-fiber 5. An analytical expression 
for the backscatter signal R bs (A<) is not available for this situation and Monte Carlo simulations 

5 were used to model the behavior of R bs (^) as a function of the diameter of the fibers 5, 6 and of 
the optical properties of the sample 1. Fig. 4 shows the results of Monte Carlo simulations using 
the MCML-code (Monte Carlo for Multi-Layered media) described in L-H Wang, S.L. Jacques, 
L-Q Zheng, "MCML - Monte Carlo Modeling of Photon Transport in Multi-Layered Tissues", 
Computer Methods and Programs in Biomedicine 47, pp. 131-146 (1995); and L-H Wang, S.L. 

10 Jacques, L-Q Zheng, "CONV - Convolution for Responses to a Finite Diameter Photon Beam 
Incident on Multi-Layered Tissues", Computer Methods and Programs in Biomedicine 54, pp. 
141-150 (1997); for a homogeneous medium with an anisotropy value g=0.9. A flat circular 
incident beam with diameter dfiber is directed onto the sample 1, and the differential backscatter 
signal R bs is calculated by subtracting the total reflectance in the c-fiber 6 (with diameter df lbe r 

1 5 and center located at a distance d flb er from the center of the incident beam) from the total 

reflectance in the dc-fiber 5 (with diameter df lbe r overlapping the incident beam). Simulations 
were performed for sets of four different scattering coefficients (|Xs=15, 25, 50 and 80 mm" 1 ), 
four different fiber diameters (d fibe r=200, 400, 600 and 800 ^im) and five different absorption 
coefficients (n a =0, 0.2, 0.4, 0.6 and 0.8 mm" 1 ). 

20 [0048] Fig. 4 shows R bs as a function of absorption coefficient where the open 

circles/dashed lines correspond to df lber =200 urn, the filled circles/dotted lines correspond to 
dfi ber =400 ^m, the open squares/solid lines correspond to dfi ber =600 urn and the filled 
squares/dash-dotted lines correspond to d fibe r=800 urn. The differential backscatter signal R bs for 
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each scattering coefficient (i s was normalized to the (dfjber = 200 |jm, ^ a =0 mm" 1 ) case. Fig. 4 
shows that in the absence of absorption, i.e. \i a = 0, the differential backscatter signal Rb S depends 
linearly on the scattering coefficient, |i s . Furthermore, the slope of the straight lines (signifying 
the relation between Rb S and n a ) depends only on the fiber diameter and is independent of the 
5 scattering coefficient \i s . The latter is more clearly demonstrated in Fig. 5, where the differential 
backscatter signal Rb S is normalized to unity at zero absorption for all scattering coefficients 
The open circles correspond to dfib e r = 200 jam and the filled squares correspond to dfiber^OO ^im. 
These Monte Carlo simulations therefore suggest that in the situation where mfp<dfiber> the 
diameter of the fibers 5, 6 determines the average path length x of the measured photons. The 

10 backscatter signal Rb S for this range of parameters can thus be written as 

K W = C\ M S - exp(-r Ju s ) = C r M s ' exp(-C 2 ■ d fiber • fu a ) (7) 
where Ci and C2 are constants, x is the average path length, |i a is the absorption coefficient, ^ s is 
the scattering coefficient and dfiber is the fiber diameter of the fibers 5, 6. 
[0049] An exact analytical expression for Rbs(k) is not available due to the large 

15 contribution of multiple scattering events to the signal. Measurements were done for 

determining a total integrated backscatter signal R to t for a range of X between 400-900 nm, using 
the formula 

900 wn 

R,M= faRn&M.) (8) 

400 nm 

[0050] Fig. 6 shows that the integrated total backscatter signal R to t(Hs) is proportional to 

20 Vl s (\) in the relevant range of, between 1 0- 1 00 mm" 1 . Therefore, in the absence of absorbers, it 
follows that 

R k sW = Ca PP -MV W 
14 



207,517 

SUBSTITUTE SPECIFICATION 

which is in agreement with the Monte Carlo simulations. 

[0051] In the presence of n absorbing substances in a suspension, with specific 

absorption coefficients Ha Spec ''0-)> the differential backscatter signal becomes 

K -C w . A (A).«p(-r. Ea'/C^W (10) 

1=1 

5 where x is the average path length of the detected backscattered photons and pi is the 
concentration of the substance i. 

[0052] Non-linear phenomena such as an inhomogeneous distribution of absorbers are 

not incorporated in Eq. (10), but can be added by the skilled person, see e.g. the van Veen et ah 
article cited herein. 

10 [0053] Fig. 7 shows the measured and calculated average path length x as a function of 

the average scattering coefficient <\i s (X)> (with 500 nm < X < 700 nm) for df lbe r = 0.4 mm and 
for absorption coefficient \i a (k) = 2.0 mm" 1 at X = 600 nm. In Fig. 7 measurements are indicated 
by dots, the x = 2 • mfp curve is indicated by a line 71 and Monte Carlo simulations are depicted 
by dashed lines. Identical results were obtained for suspensions with an absorption coefficient of 

1 5 Ha=l 0 mm" 1 at 600 nm. The average path length x was determined using suspensions of 

polystyrene spheres with different sizes and concentrations to vary the scattering coefficient \i s (k). 
The anisotropy g of these suspensions was in the range of 0.8-0.9. Evans Blue dye was added as 
an absorber, and the average path length x was calculated from Eqs. (9) and (10) and knowledge 
of the concentrations and a specific absorption coefficient of Evans Blue, as will be known to the 

20 skilled person. 
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[0054] Looking at the measured average path lengths of Fig. 7, it clearly shows that for 

large scattering coefficients (|i s =10-100 mm" 1 the range relevant for tissue) the average path 
length x is independent of the scattering coefficient \i s to within 10% and approximately equal to 
half the fiber diameter (x ~ 0.24 mm while d fib er = 0.40 mm), in agreement with the Monte Carlo 

5 simulations (the dashed lines correspond to Monte Carlo calculations for dfiber = 0.2, 0.4, 0.6 and 
0.8 mm). For small scattering coefficients (e.g. \i s < 5 mm" 1 ) the average path length x is well 
described by x = 2 • mfp, according to Eq. (4), see line 71. Fig. 7 also clearly demonstrates that 
the transition from the 'single scattering regime' to the 'constant path length regime' occurs for 
mean free paths of the order of the fiber diameter. It is therefore expected that single scattering 

10 will prevail over a larger range of scattering coefficients for fiber diameters smaller than 400 \xm. 
[0055] In the following, the effect of absorption on the average path length x will been 

examined in more detail. Various concentrations of Evans Blue dye were added to a suspension 
of polystyrene spheres with scattering coefficients fa s of 35 mm" 1 . The concentrations of Evans 
Blue (EB) dye were varied such that the absorption coefficient jx s at 600 nm was in the range of 0 

15 to 2 mm' 1 . Typical results of the differential backscatter signal R bs for three different absorption 
coefficients \i a are shown in Fig. 8. Note that the signal I of the dc-fiber 5 is plotted on a 
different vertical scale than the signal J of the c-fiber 6 and the differential backscatter signal Rb S . 
[0056] The spectra with Evans Blue R EB present in the suspension were divided by the 

spectrum with no Evans Blue present R° and the negative natural logarithm of the ratio R EB / R° 



20 was determined: 



A = -ln(R EB /R°) = T-p-Mr' EB 00 



where p is the concentration of Evans Blue, and ^i a spec ' EB is the specific absorption coefficient of 



Evans Blue. 
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[0057] Fig. 9 shows a typical spectrum of an absorption curve 92, along with the specific 

absorption coefficient ja a spec,HB of Evans Blue dye, see curve 94. 

[0058] For all concentrations, an area A* under the absorption curve 92 was determined 

in the wavelength range X between 500 and 650 nm. From Eq. (1 1) it follows that if the average 
5 path length x is independent of the absorption coefficient |a a spec > area A* should depend linearly 
on the concentration p of the species in the suspension. 

[0059] Fig. 10 shows a measured A* as a function of the absorption coefficient \x a at 600 

nm. A curve fitting is done, resulting in a line 104 for ii s (X) = 35 mm' 1 . Fig. 10 shows that the 
average path length x is indeed independent of the absorption coefficient (Li a in the range 0-2 
10 mm" 1 . 

[0060] From the previous results of Figs. 4 to 10, it shows that for mfp<d fi ber, the 

differential backscatter signal Rb S is described by Eq. (7) with C 2 ~ 0.6. Fig. 1 1 shows typical 
spectra measured in a suspension of 1.0 (im polystyrene spheres with and without Evans Blue 
dye (^i a = 2 and 0 mm" 1 at 600 nm, respectively). From Eq. (7) and Fig. 7, the relation between 
1 5 the differential backscatter signals with and without absorber is given by 

R bs (A. Ma ) = R ts W>) ■ exp(-0 .24 • Ma ) (12) 
[0061] The calculated spectrum according to Eq. (12) is plotted as a dashed line 1 10 in 

Fig. 1 1 and shows excellent agreement with the measured R bs (k,Ha), see line 1 1 1. In Fig. 11, line 
112 depicts R bs (^, 0). 

20 [0062] In short, the average path length x of photons measured when subtracting the 

signals of the c-fiber 6 from the dc-fiber 5 using Eq. (1) is independent of the optical properties 
of the sample 1 and approximately equal to half the diameter of the fibers 5, 6 used, as long as 
the fiber diameter is larger than the mfp. 
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[0063] In a specific embodiment, the device according to the invention is arranged to 

determine concentrations of oxygenated blood in tissue. Since the scattering coefficient of tissue 
^ s tissue is in the range of 10-100 mm" 1 , the fiber diameter should be smaller than a certain 
maximum diameter d max where d max is between 10 and 100 jim, for example, smaller than 50 nm, 
5 in order to measure predominantly single scattering in tissue. In this case, Eq. (5) holds. For 
fibers 5, 6 with much larger diameters (e.g. 200 or 400 ^im), the differential backscatter signal 
R bs (k) is described by Eq. (10) with x ~ 0.6 ■ df,b er . 

[0064] It is presently known that the wavelength dependence of the scattering coefficient 

in tissue ^i s tissue can be adequately described by an empirical power-law function, see also the van 
10 de Hulst article cited herein, as well as A. Angstrom, "On the Atmospheric Transmission of Sun 
Radiation and on Dust in the Air", Geograf. Ann. Deut. 1 1, pp. 156-166 (1929); and R. Graaff, 
J.G. Aarnoudse, J.R. Zijp, P.M. A. Sloot, F.F.M. de Mul, J. Greve, and M.H. Koelink, "Reduced 
Light-Scattering Properties for Mixtures of Spherical Particles: A Simple Approximation 
Derived from Mie Calculations", Applied Optics. 31, pp. 1370-1376 (1992). 

15 mT u€ W = a ' X b < 13 ) 

with a and b constants that depend on the size, concentration and relative refractive index of the 
scatterers (i.e. substances) present in the detection volume. 

[0065] The dominant absorbers in tissue in the visible wavelength range are oxygenated 

and deoxygenated blood. Thus in tissue Eq. (10) becomes 

R b M) = C m ■ aX-» ■ exp(-0.6 • d fiber • p blood ■ (S 02 ■ fi-* + d - S 02 ) • (14) 
= C W -X b .exp(-0.6.^ r - Pblood <S 02 - M r ox +Q-S 02 ) »r J ' aa )) 

where p b | 00 d is the concentration of blood, S02 is the blood oxygenation (percentage oxygen 
saturation) in a certain detection volume, C app is a constant that depends on the calibration 
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constant c, C' a pp is C app -a, X is the wavelength, b is the slope of the scattering coefficient defined 
in Eq. 13, n a spec,ox is the specific absorption coefficients of fully oxygenated blood, and ^ s P ec ' deox 
is the specific absorption coefficients of fully deoxygenated blood. 

[0066] Non-linear phenomena such as an inhomogeneous distribution of absorbers are 

not incorporated in Eq. (14), but can be added by the skilled person, 

[0067] Since the specific absorption coefficients of fully oxygenated (^ a spec '° x ) and fully 

deoxygenated (|i a spec ' deox ) blood are well known, see Fig. 12, Eq. (14) can be fitted to the 



measured data to yield the slope b of the scattering coefficient ^ s 18 ue , the concentration pbiood and 
the oxygen saturation S02 of the blood present in the detection volume. When a correction is 

10 made for the inhomogeneous distribution of blood in the vessels, a vessel diameter D may be 
determined as well. Since the average detection depth is small (e.g. 0.1 mm), the blood present 
in the detection volume when measuring non-invasively is located in capillaries. 
[0068] In Fig. 13, in vivo measurements of backscattering in a human trachea together 

with a fit using Eq. (14) are shown. The measurements are realized using a fiber diameter of 400 

15 |im. The dots depict the measurements and curve 130 is a fitting curve. In Fig. 13, b =-0.94 and 
the oxygenation So2 = 95%. 

[0069] The present invention can be used for tumor detection. Tumor growth may, due 

to its excessive oxygen consumption, be accompanied by a low capillary oxygen saturation, 
which can only be assessed using a very localized measurement. Since (pre-)cancerous tissue is 
20 generally more heterogeneous than normal tissue, the standard deviation of multiple 

measurements is likely to be larger for (pre-)cancerous tissue than for normal tissue. Standard 
deviations in the measurements can be calculated for the oxygen saturation, the blood 

tissue 

concentration, the blood vessel diameter, and the slope b of the scattering coefficient n s issue . It is 
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noted that the invention is by no means restricted to determine a concentration of a substance as 
the physical feature. All features, mentioned in the previous phrase can be regarded as physical 
features. 

[0070] An example of a measurement of a lung tumor is shown in Fig. 14. The shape of 

5 the dip in the wavelength range of 500-600 nm in this figure demonstrates the depletion of 

oxygen from the capillaries of this tumor due to its excessive oxygen consumption. 

[0071] When a needle-probe is used, the local oxygenation and scattering coefficient 

Hs tlssue , can be measured invasively. This could be helpful in determining tumor-margins intra- 

operatively in real-time, for instance during resection of a breast-tumor. 
10 [0072] According to an embodiment, the device comprises multiple probes and a 

multichannel spectrometer for multiple simultaneous measurements on different locations of the 

sample 1 . Using this device, multiple measurements can be made simultaneously on different 

locations of for example a suspicious lesion. 

[0073] In yet another embodiment, the device comprises at least two pairs of fibers, 

1 5 having different fiber diameters. For example, when a pair of fibers with 1 00 jim, a pair of fibers 
with a diameter of 200 ^im and a pair of fibers with a diameter of 400 |im are used, information 
from different depths in the sample 1 can be obtained as the average path length increases with 
increasing fiber diameter. 

[0074] The method and apparatus according to the invention can also be used to analyze 

20 local drug concentrations. From Eq. (10) it follows that if the specific absorption coefficient of a 
certain drug is known, the local concentration p of that substance can be determined using the 
invention. 
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[0075] Another possibility of the present invention is to monitor glucose concentrations. 

The scattering coefficient n s tissue depends among others on the relative refractive index of the 
scatterers with respect to the surrounding medium (in tissue: cytoplasm). The refractive index 
of the surrounding cytoplasm is likely to depend on the concentration of glucose. A change in 
5 the glucose concentration will therefore likely affect the slope b of the scattering coefficient 
^i s tissue , see Eq.(13). 

[0076] While specific embodiments of the invention have been described above, it will 

be appreciated that the invention may be practiced otherwise than as described. For example, a 
concentration of a substance in polluted water may be calculated. The description is not intended 
10 to limit the scope of the invention. 
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